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Abstract 

A  three-degree-of-freedom  mathematical  model  has  been  developed  to  describe  the  zigzag 
manoeuvring  behaviour  of  a  prototype  hull  form  used  in  the  Royal  Australian  Navy 
precursor  mine  sweeping  drone  boat.  The  model  predicts  the  yaw  rale  of  the  vessel  in 
response  to  the  angle  of  thrust  that  a  single  outboard  motor  makes  urith  respect  to  the 
centre  line  of  the  boat.  The  outboard  motor  is  assumed  to  supply  a  constant  thrust.  This 
report  describes  the  development  of  the  model  from  first  principles  and  its  transformation 
into  a  form  suitable  for  use  in  closed  loop  control  modelling.  The  coefficients  of  the  model 
were  determined  using  full  scale  manoeuvring  data  for  the  prototype  hull  form. 
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Development  of  a  Manoeuvring  Model 
for  the  RAN  Precursor  Mine  Sweeping 

Drone  Boat 


1.  Introduction 

The  Precursor  Mine  Sweeping  Drone  Boat  under  development  by  the  Royal 
Australian  Navy  will  be  based  on  a  commercially  available  tri-hulled  vessel  made 
of  fibreglass.  The  vessel  will  be  eight  metres  in  length  and  powered  by  two  1000 
kW  outboard  motors.  The  boat  will  be  fitted  with  a  radio  position  fixing  system 
and  will  tow  an  acoustic  influence  mine  sweeping  device. 

The  mine  sweeping  drone  boat  will  be  operated  remotely  via  a  radio  link  to  a 
shore  station.  The  boat  will  be  equipped  with  a  programmable  autopilot  to 
maintain  the  vessel  dose  to  a  defined  track.  The  drone  boat  concept  is  depicted  in 
figure  (1). 


Figure  1;  The  drone  boat  towing  a  mine  sweeping  device.  (The  sweeping  device  is 
submerged,  the  semi-submerged  towfish  indicates  its  position.) 


The  work  described  in  this  report  was  conducted  at  a  time  prior  to  the 
construction  of  the  drone  boat  hull.  However,  the  RAN  supplied  a  six  metre 
version  of  the  same  hull  form  which  was  used  for  preliminary  investigations  of 
the  concept.  This  report  describes  the  development  of  a  horizontal  plane,  three- 


degree-of-freedom  mathematical  manoeuvring  model  for  this  prototype  of  the 
drone  boat  hull.  The  model  predicts  the  yaw  rate  of  the  vessel  as  a  function  of  the 
angle  of  thrust.  The  outboard  motor  is  assumed  to  supply  a  constant  thrust.  The 
effects  of  wind  and  sea-state  are  not  considered.  The  manoeuvring  model  was 
specifically  developed  to  simulate  and  compare  the  performance  of  two 
commercially  available  autopilot  systems  (Campanella,  Anderson  and  Spirovski). 

The  specification  for  the  drone  boat  autopilot  system  is  shown  in  figure  2.  The 
drone  boat  autopilot  is  programmed  with  baseline  coordinates  defined  by  points 
A  and  B.  The  distance  the  boat  wanders  from  this  baseline  is  called  the  cross-track 
error  ( Yte ).  This  must  not  be  greater  than  ten  metres  on  either  side  of  the  baseline. 
The  autopilot  control  law  is  based  on  measurements  of  the  instantaneous  heading 
of  the  boat  (y)  with  respect  to  the  track,  the  known  required  track  heading  or 
reference  heading,  and  the  cross-track  error  (Yu). 


Maximum  cross  track  error 


Figure  2:  Guidance  and  control  specification  for  the  drone  boat  autopilot  system. 


2.  Derivation  of  the  Manoeuvring  Model 

The  manoeuvring  model  is  required  to  predict  the  vessel  rate  of  change  of 
direction  (yaw  rate,  r)  in  response  to  the  change  in  angle  of  the  outboard  motor 
with  respect  to  the  centreline  of  the  boat  ( rudder  angle,  § ).  The  outboard  motor 
is  assumed  to  produce  constant  thrust.  The  model  is  developed  for  calm  water 
and  neglects  forces  on  the  boat  due  to  wind  and  current.  The  tow  device  has  been 
assumed  to  be  a  rigid  part  of  the  drone  boat.  The  extra  drag  force  due  to  the  tow 
is  accounted  for  in  the  model  coefficients. 


Figure  3:  A  body  fixed  coordinate  system  for  surface  vessels. 

Figure  3  shows  the  coordinate  system  used  in  the  equations  of  motion  for  the 
drone  boat,  where  X0  and  yo  are  the  global  coordinate  axes  fixed  in  space  and 
where  xo  lies  in  the  same  direction  as  the  baseline  AB.  The  forward  direction  of 
the  boat  lies  along  the  x-axis  while  the  y-axis  denotes  sidewards  direction.  The  x 
and  y  axes  are  centred  at  a  fixed  arbitrary  point  O  on  the  boat  to  form  the  local  or 
body  fixed  coordinate  system.  This  coordinate  system  follows  the  conventional 
right  hand  rule,  giving  all  rotations  a  positive  sense  in  the  clockwise  direction. 
The  angle  \|/  represents  the  angle  of  rotation  of  the  boat's  x-axis  with  respect  to 
the  global  x0-axis. 
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HEAVE 


SURGE 


Figure  4:  Displacement  motions  of  a  surface  vessel. 


Figure  4  illustrates  the  three  displacement  motions;  heave,  sway  and  surge.  The 
motion  of  the  boat  in  the  z  axial  direction  is  designated  heave  and  is  positive 
downwards.  Sidewards  motion  of  the  boat  is  sway.  Surge  is  the  third  linear 
motion  and  is  due  to  forward  displacement  in  the  x -axial  direction.  Heave  is 
neglected  in  the  horizontal  manoeuvring  model. 

Figure  5  illustrates  the  three  rotational  motions;  pitch,  roll,  and  yaw. 

Roll  is  the  rotation  about  the  x-axis,  while  rotation  about  the  y-axis  is  pitch.  Yaw 
is  the  rotation  about  the  z-axis  of  the  vessel.  Pitch  and  roll  are  neglected  in  the 
following  horizontal  plane  model.  In  theory,  any  of  the  six  motions  can  occur  in 
the  absence  of  the  other,  but  in  practice  coupled  motions  will  exist. 
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Figitre  5:  Rotational  motions  of  a  surface  vessel. 


2.1  Mathematical  Development 

The  following  is  a  derivation  of  the  manoeuvring  equations  for  the  mathematical 
model,  based  on  the  work  of  Abkowitz  (1969),  Gill  (1979)  and  Clarke  (1982).  The 
essential  assumptions  of  the  model  are  summarised  as  follows: 

1.  The  variation  in  the  mass  of  the  drone  boat  due  to  the  consumption 
of  fuel  is  negligible. 

2.  Forces  on  the  boat  due  to  the  rotation  of  the  Earth  are  neglected. 

This  is  equi  valent  to  neglecting  Coriolis  accelerations. 

3.  The  model  is  linear. 

4.  Pitch,  roll,  and  heave  are  neglected. 

5.  Speed  loss  in  the  surge  direction  is  negligible. 

6.  The  thrust  produced  by  the  propellor  is  constant  and  moments  due 
to  the  rotation  of  the  propellor  are  neglected. 

7.  The  boat  is  symmetrical  about  the  x*a xis. 

8.  The  boat/tow  system  is  assumed  rigid  in  the  model. 

9.  The  oscillatory  nature  of  the  manoeuvres  is  assumed  to  have 
negligible  effect  on  the  external  forces  applied  to  the  boat/  tow 
system. 

Consider  the  vessel  in  figure  (3)  proceeding  along  some  arbitrary  path  with  linear 
and  rotational  velocity  UG  and  flG  respectively.  The  application  of  Newton's 
second  law  of  motion  expressed  in  the  global  coordinate  system  yields  the 
following  equations: 
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F“  j-(mUG) 

dt 

...a) 

dt 

...(2) 

Where  F  =  The  resultant  externally  applied  force  vector  acting  on  the  boat 
at  the  centre  of  gravity  (considered  as  a  column  vector). 

Mg  3  The  resultant  externally  applied  moment  vector  acting  ©?»  the 

boat  about  the  centre  of  gravity  (considered  as  a  column  vector). 

UG  *  The  linear  velocity  vector  with  respect  to  the  centre  of  gravity 
(considered  as  a  column  vector). 

Slc  =*  The  rotational  velocity  vector  of  the  boat  about  the  centre  of 
gravity  (considered  as  a  column  vector). 

m  3  The  dry  mass  of  the  body. 

Ig  “  A  3x3  matrix  of  the  body's  moment  of  inertia  about  an 

orthogonal  axes  parallel  to  x,y,z  through  the  centre  of  gravity. 

Although  the  expressions  are  more  complicated,  modelling  the  forces  about  the 
local  origin  O  is  desirable  since  O  can  be  placed  at  a  point  on  the  body  where 
symmetry  might  be  considered  to  exist.  Figure  (6)  shows  the  vector  relationship 
between  the  two  coordinate  systems,  one  located  at  the  centre  of  gravity,  and  the 
other  at  an  origin  O. 


u= 


Figure  6:  Two  body  fixed  coordinate  systems  showing  the  body's  motion  vectors. 

Since  the  velocity  at  the  centre  of  gravity  UG  is  equal  to  the  velocity  U  at  the 
origin  O,  plus  the  velocity  of  the  centre  of  gravity  relative  to  O  i.e. 

UG  =  U  +  QxR,  then  equation  (1)  can  be  expressed  as  shown  in  equation  (3). 
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F  =  m— (U+QxR) 
dt 


...  (3) 


Where  R  = 

Cl 

U 

x  = 


The  position  vector  of  the  boat  centre  of  gravity  with 
respect  to  the  origin  O. 

The  rotational  velocity  vector  of  the  boat  about  the  origin  O. 
The  linear  velocity  of  the  boat  at  the  origin  O. 

The  vector  cross  product  operator. 


The  moment  about  the  origin  O  is  equal  to  the  moment  about  the  centre  of  gravity 
plus  the  moment  due  to  the  force  F  acting  over  the  distance  R.  The  moment  is 
therefore  given  by 

M  =  Mg  +  RxF 


Substituting  for  and  F,  where 


..  d  ,  .  (  dV 

dt  dt2  J 

and  F  is  given  by  equation  (3)  the  following  expression  for  the  moment  about  the  origin 
results. 


M«— (Kl)+mRx—  (U) 
dt  dt 

where  I  is  a  3x3  matrix  of  the  moment  of  inertia  about  the  origin  O, 


...(4) 


such  that  I a 


Ixx 

_I*y 

-Ixz 

^yx 

_Iyz 

Izx 

'lzy 

Izz 

and 


M  is  the  moment  acting  on  the  boat  about  the  centre  of  gravity 
(considered  as  a  column  vector). 

Since  pitch,  heave  and  roll  are  neglected,  the  vectors  arising  in  the  model  are: 


F 

= 

(X,  Y,  0)T 

M 

= 

(0,  0,  n)t 

« 

s 

(0,  0,  r)T 

U 

« 

(u,  v,  0)T 

R 

= 

(xG,  0,  0)T 
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Where  u  =  The  component  of  the  boat's  velocity  resolved  along  the  x-axis. 

v  =  The  component  of  the  boat's  velocity  resolved  along  the  y-axis. 

r  =  The  component  of  angular  velocity  of  the  boat  resolved  about 

the  z-axis. 

x0  =  The  position  of  the  boat's  centre  of  gravity  along  the  x-axis 

with  respect  to  O. 

Expanding  equations  (3)  and  (4)  into  their  components,  and  using  the  values  given 
for  F,  M,  Q,  U  and  R  above,  results  in  the  following  expressions: 

m(u-rv-xGr2)  =  X  ...(5) 

•  • 

m(v+ru  +  xG  r)  =  Y  ...(6) 

•  • 

1 22  r  +  mxG(v+ru)  =  N  ...(7) 

Now  consider  the  vessel  to  have  velocities  and  accelerations  given  by 
uo've'ro'Uo,Vo,r0If  each  term  is  perturbed  by  a  small  amount  then  the 
resultant  velocities  and  accelerations  are  given  by: 

u  =  Au  +  u0 
v  =  Av  + vQ 
r  =  Ar+r0 

•  •  • 

u  =  Au+u0 

•  mm 

V  =  A  v+  v0 

•  •  • 

r  =  Ar+r0 

Y  =  Y0  +  AY 


Consider  now  the  development  of  the  sway  force  model.  Expanding  the  sway 
force  equation  (6)  in  terms  of  the  above  variables  gives: 

ml(Av+v0)+(Ar  +  r0)(Au  +  u0)  +  xG(Ar+r0)]  =  Y0+AY  ..(8) 

Equation  (8)  is  simplified  by  considering  the  case  of  steady  state  motion  where  the 
vessel  is  travelling  with  constant  forward  speed,  with  small  disturbances 
occurring  about  this  motion.  Hence 

•  •  • 

v0  =  v0  =  rD  =  r0  =  u0  =  Y0  =  0,  and  U,,  *■  0. 

•  •  •  I  •  • 

Therefore  v  =  Av,r  =  Ar,u  =  A  u,  v  =  A  v,r  =  Ar,Y  =  AY  and  u  =  Au  +  uc. 
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Neglecting  the  non-linear  terms,  the  expression  for  the  sway  force  is  given  by; 


m(v+ruc+xGr)-Y  ...(9) 

Similarly,  the  yaw  moment  equation  is  given  by 

Izz  r+mxG(v+ru0)*N  ...(10) 

and  the  surge  equation  is  given  by 


m(u)«X  ...(11) 

Assuming  that  Au  is  small,  then  U  is  negligible,  therefore  equation  (11)  is 
neglected  here.  Hence  equations  (9)  and  (10)  form  the  basis  of  the  mathematical 
model. 

X  and  Y  denote  the  sum  of  resolved  components  of  external  forces  in  the  surge 
and  sway  directions  respectively,  and  N  gives  the  sum  of  external  moments  about 
the  z-axis.  In  the  case  of  the  net  external  sway  force: 

v-Iy, 

=  yt+yh+ye 

Where  YT  are  the  thrust  forces  in  the  sway  direction, 

Yh  are  the  hydrodynamic  forces  in  the  sway  direction,  and 
Ye  are  any  other  external  forces  which  act  on  the  body. 

Assuming  that  YE  =  0,  the  hydrodynamic  and  thrust  forces  which  act  on  the 
drone  boat  in  the  sway  direction  are  assumed  to  be  a  function  of  the  variables 
v,  v,r,r,8,  such  that: 


y  =  y„+yt 

=  Y(v,v,r,r,8 


...(12) 


Where  5  is  the  angle  of  the  control  surface  (or  in  the  case  of  the  drone  boat,  the 
outboard  motor)  with  respect  to  the  centreline  of  the  body.  5  is  the  resultant  of  a 
steady  value  50and  a  small  increment  A5. 

A  first  order  Taylor  series  expansion  of  the  sway  force  is  formed  about  the 
equilibrium  values  of  the  dynamic  variables  and  neglecting  terms  higher  than  the 
first  order,  giving 


3Y  3Y  •  3Y  9Y  •  3Y 

Y  “  Y  + —  Av  +  —  A  v+ —  Ar  + —  A  r+ —  AS 


...  (13) 


d  v, 


3vC  ' 


3r, 


drWte 


98  tmt0 


The  literature  expresses  the  partial  derivatives  as  follows: 


Y 

v  3vv 


etc. 


This  is  not  to  be  confused  with  the  variables  YT  and  YH  which  use  uppercase 
letters  for  the  subscripts  to  denote  a  particular  component  of  the  sway  force. 
Using  this  notation  and  recognising  that  the  steady  state  sway  force  Y0  =  0,  and 
the  steady  state  outboard  motor  angle  S0*  0,  equation  (13)  becomes 


Y  =  Yvv  +  Y.  v+Yrr+Y.  r+Y,8 

v  r 


...  (14) 


The  term  YsS  models  the  sway  force  of  the  outboard  motor's  thrust.  This  is  a  very 
simple  model  which  assumes  that  the  motor  produces  constant  thrust  and  that  5 
is  small. 

The  partial  derivatives,  otherwise  known  as  hydrodynamic  coefficients  or 
derivatives,  are  a  function  of  frequency.  Therefore,  the  correct  representation  for 
Yv  is  Yv  (to)  where  (0  is  the  angular  frequency  associated  with  the  oscillation  of 
the  boat's  motion.  The  frequency  of  the  manoeuvres  for  the  drone  boat  is  expected 
to  be  small  and  so  the  coefficients  can  be  taken  as  constants.  Substituting  equation 
(14)  into  (9)  and  rearranging,  gives: 


(Y.  -m)v+Yvv+(Y.  -mxG)r+(Y, -mu0)  r+Ys8  =  0  ...(15) 

v  r 


The  same  procedure  applied  to  the  yaw  moment  equation  (10)  gives: 

(N.  •mxG)v+Nlv+(N,  -  Ilz)r+(N, -mxG uo)r  +  Ns8  =  0  ...(16) 

v  r 


The  manoeuvring  model,  defined  by  equations  (15)  and  (16)  was  used  as  the  open 
loop  transfer  function  to  be  controlled  by  the  drone  boat  autopilot  control  law.  For 
ease  of  computer  simulation  the  open  loop  manoeuvring  model  and  closed  loop 
control  model  were  transformed  to  the  z-domain.  An  s-domain  representation  is 
required  as  an  intermediate  step  (Campanella,  Anderson,  and  Spirovski). 

The  differential  operator  d/dt  is  replaced  with  the  Laplacian  s,  in  equations  (15) 
and  (16),  and  rearranged  to  give  the  following  equation: 

I  (Es  +  F)I  -(As  +  B)J 

8  (As  +  B)(Cs  +  D)-(Es+F)(Gs  +  H) 

where  r  =  r(s)  [the  Laplace  transform  of  r(t)J 

8  =  8(s)  [the  Laplace  transform  of  8(t)J 

A  =  Y.  -m 


B  =  Yv 
C  =  N.-I„ 

r 
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D  »  Nr  -  mxc  uB 
E  ■  N.  -mxc 


F  ■  Nv 
G  -  Y.  -mxG 

r 

H  -  Y,  -  mu0 
I  =  Y* 

J  -N, 

Equation  (17)  can  be  simplified  further  to  give: 

I .  ...os 

8  sz+bs+c 

Where  5  and  r  are  the  outboard  angle  and  the  yaw  rate  in  the  s-plane 
respectively.  K,  a,  b,  and  c  are  functions  of  the  hydrodynamic  coefficients.  This 
equation  represents  the  open  loop  transfer  function  given  in  the  s-domain  of  the 
change  in  boat  heading  in  response  to  change  in  angle  of  thrust  of  the  outboard 
motor. 


3.  Determining  the  Model  Coefficients 

The  model  parameters  K,  a,  b,  and  c  were  determined  by  fitting  system  response 
measurements  to  the  model  form  described  in  equation  (18).  The  process  of 
developing  the  mathematical  manoeuvring  model  of  a  particular  vehicle  is 
summarised  in  figure  (7). 

The  formulation  of  the  model  can  be  approached  in  either  of  two  ways;  firstly 
using  a  physical  approach  which  builds  a  model  based  on  known  properties  and 
an  understanding  of  the  underlying  physical  processes  of  the  system,  or  secondly 
using  a  non-physical  approach  where  the  physical  properties  or  mechanisms  of 
the  system  may  not  be  completely  known,  or  are  too  complex  to  be  modelled,  or 
both. 

The  approach  used  in  this  report  was  a  combination  of  both  in  that  the  form  of 
the  model  was  derived  from  force  and  moment  equations  expressed  with 
hydrodynamic  coefficients  in  a  body-fixed  coordinate  system.  For  ease  of 
computer  simulation  the  expressions  were  transformed  to  a  transfer  function  in 
the  s-domain  where  the  coefficients  do  not  relate  simply  to  the  standard 
hydrodynamic  derivatives. 

There  are  several  approaches  to  the  evaluation  of  hydrodynamic  coefficients. 
One  such  approach  is  the  tank  testing  of  scale  models.  Scale  model  testing 
involves  subjecting  a  scale  model  of  the  hull  form  to  precisely  controlled  motions 
in  a  fluid  flow.  Forces  acting  on  the  scale  model  are  measured  and  from  this 
hydrodynamic  coefficients  are  derived.  This  method  involves  the  use  of 
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specialised  facilities  such  as  towing  tanks,  circulating  water  channels,  and  planar 
motion  mechanisms  (Gill  1979). 

The  approach  used  here  was  to  derive  the  required  modelling  parameters  from 
full-scale  manoeuvres.  In  its  most  basic  form,  this  involves  conducting 
manoeuvres  with  the  actual  vessel,  in  still  waters,  while  measuring  control  surface 
angles  (the  stimulus),  and  the  boat  heading  (the  response).  Recursive  algorithms 
or  trial  and  error  methods  can  be  employed  to  fit  the  experimental  data  to  the 
postulated  model  form  and  so  derive  the  model  parameters. 
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Figure  7:  The  system  identification  process. 
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4.  Full-Scale  Drone  Boat  Manoeuvring 
Experiment 

A  manoeuvring  trial  of  the  prototype  drone  boat  was  conducted  at  HMAS 
Crowell,  Jervis  Bay,  Australian  Capital  Territory.  Set  manoeuvres  were  performed 
and  measurements  of  rudder  angle,  boat  heading,  and  speed  in  the  surge 
direction  were  taken.  Manoeuvres  were  conducted  with  and  without  the  towed 
acoustic  mine  sweeping  device.  All  manoeuvres  were  performed  in  zero  sea-state 
by  a  helmsman. 

Renilson  (1968)  and  Abkowitz  (1969)  describe  several  types  of  manoeuvres 
suitable  for  characterising  the  full  scale  manoeuvring  performance  of  vessels.  The 
"zigzag*  manoeuvre  was  used  in  this  experiment.  As  a  transient  manoeuvre  it 
allows  both  the  acceleration  and  velocity  components  to  be  derived.  Abkowitz 
(1960)  considers  a  10°  zigzag  to  be  a  mild  manoeuvre,  where  the  significant 
contribution  to  the  forces  is  made  by  linear  terms,  which  is  appropriate  for  the 
drone  boat  manoeuvring  model. 

In  order  to  perform  the  zigzag,  the  boat  was  set  on  a  straight  course  at  a 
constant  initial  velocity  of  7  knots.  For  a  10s  manoeuvre,  the  helmsman  set  the 
rudder  10°  to  starboard.  When  the  boat  heading  changed  by  10°,  the  rudder  was 
then  set  10°  to  port  By  continuing  this  action  the  vessel  was  made  to  zigzag. 

Manoeuvring  data  were  recorded  for  both  a  10°  and  a  5°  zigzag.  The  acoustic 
mine  sweeping  device  was  attached  for  a  10*  manoeuvre. 

The  boat  was  fitted  with  a  Coursemaster  Model  CM500  Series  2  marine 
autopilot  system  which  was  primarily  used  as  the  means  of  digitising  and  scaling 
the  boat  heading  sensor  signal.  The  rudder  angle  and  boat  heading  sensors  were 
manufactured  and  fitted  by  Coursemaster  for  the  support  of  their  autopilot 
system.  Boat  speed  was  measured  with  a  conventional  impellor-speed-log. 

Rudder  angle  was  measured  with  a  potentiometer-based  transducer.  A  fluxgate 
compass  was  used  for  measuring  the  boat  heading. 

The  digitised  boat  heading  data  were  transferred  to  an  MRL-developed  data 
recorder.  The  drone  boat  was  fitted  with  a  HamilHaven  Type  8904  data  logger 
and  this  was  used  to  acquire  velocity  and  rudder  angle  data.  Data  were 
transferred  to  the  MRL  data  logger  after  each  manoeuvre. 

The  MRL  data  recorder  was  a  Toshiba  T1 600/40  IBM  compatible  "laptop" 
personal  computer.  Data  acquisition  software  was  developed  in  Turbo  Pascal 
which  read  data  from  the  computer  serial  port  and  wrote  them  to  a  disk  file. 


5.  Results 

A  graphical  computer  simulation  of  the  manoeuvring  model  was  developed  in 
DOS  Turbo  Pascal  and  implemented  on  an  IBM  compatible  personal  computer. 
The  input  to  the  program  requires,  rudder  angle  as  a  function  of  time,  the  boat's 
initial  forward  speed,  and  a  set  of  values  for  the  coefficients  described  in  equation 
(18).  The  program  predicts  the  boat's  response  to  the  input  data  and  plots  boat 
heading  as  a  function  of  time.  The  values  of  the  coefficients  shown  below  were 
derived  by  trial  and  error.  The  comparison  between  the  simulation  and  the  trial 
data  allowed  the  coefficients  to  be  adjusted  until  the  best  fit,  as  judged  by  eye 
were  obtained.  Initial  values  for  the  coefficients  were  generated  from  estimates  of 
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Boat  Heading,  rad 


hydrodynamic  derivatives  for  the  drone  boat  hull  form  supplied  by 
Renilson(1990). 

Figures  (8)  and  (9)  show  simulated  and  actual  boat  heading  for  two  different 
manoeuvres  of  the  drone  boat  without  a  tow.  Both  manoeuvres  were  performed 
with  an  initial  forward  speed  of  7  knots. 
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Figure  8:  A  5°  zigzag  manoeuvre. 

The  model  coefficients  used  to  generate  the  simulated  boat  headings  in  figures  (8) 
and  (9)  were: 

K  -  -0.984 

a  -  0.109 

b  -  0.159 

c  -  1.692 

Figure  (10)  shows  the  comparison  for  a  10°  zigzag  manoeuvre  for  the  drone  boat 
with  its  tow  load  and  an  initial  forward  speed  of  seven  knots. 

The  coefficients  for  this  simulation  were  : 

K  -  -0.555 

a  -  0.110 

b  -  0.141 

c  -  1.120 


Boat  Heading,  rad  Boat  Heading,  rad 


Time,  s 


Figure  9r  A  10°  zigzag  manoeuvre. 


Figure  10:  A  VF  zigzag  manoeuvre  performed  while  towing  the  pipe  noise  maker. 
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6.  Discussion 


The  literature  describes  various  forms  for  mathematical  manoeuvring  models.  The 
form  used  in  this  paper  (equation  17)  as  described  by  Clarke(1982),  was  found  to 
be  the  most  suitable  for  its  application  to  the  modelling  of  dosed  loop  control 
laws. 

The  zigzag  is  expected  to  be  the  predominant  manoeuvre  the  drone  boat  will 
perform  whilst  recovering  from  a  disturbance  under  the  control  of  its  autopilot. 
Hence  the  zigzag  manoeuvre  was  deemed  to  be  the  most  appropriate  motion  to 
validate  the  model. 

Figures  (8)  and  (9)  show  dose  agreement  between  simulation  and  trial  data  for 
two  different  zigzag  manoeuvres  with  no  tow  load  attached  to  the  boat.  The 
assumption  that  the  boat  and  towed  load  can  be  described  by  the  same  model 
form  appeared  to  be  justified  by  the  result  obtained  in  figure  (10). 

At  the  time  this  work  was  conducted  the  parameters  were  required  urgently  for 
the  subsequent  simulation  of  autopilot  proposals  for  the  RAN  and  so  for  speed  of 
response  the  model  parameters  were  determined  by  trial  and  error.  Whilst  this 
was  a  somewhat  tedious  process,  it  did  allow  for  the  provision  of  timely  advice. 
Trial  and  error  appeared  to  be  a  satisfactory  process  for  the  simple  model 
developed.  However,  the  procedure  can  be  automated  by  implementing  a  least 
squares  method  or  by  more  advanced  recursive  methods  described  by 
Ljung(1986).  This  would  be  especially  necessary  for  more  complex  models  with  a 
greater  number  of  coefficients. 

In  all  three  comparisons  it  can  be  observed  that  the  greatest  discrepancy 
between  simulation  and  actual  data  occurred  at  the  points  where  the  boat  changes 
its  direction.  This  error  can  be  attributed  to  the  fact  that  the  model  was  developed 
from  linear  equations  which  assume  that  perturbations  from  the  equilibrium  state 
are  small.  Under  this  assumption  the  surge  equation  is  neglected.  Closer 
agreement  might  be  obtained  by  including  the  surge  equation  in  the  model. 

Notwithstanding  the  encouraging  results  obtained  it  should  be  emphasized  that 
the  model  was  formulated  to  describe  a  boat  travelling  at  some  forward  constant 
speed  whilst  zigzagging  due  to  small  changes  in  the  angle  of  thrust  from  a  single 
outboard  motor.  This  simple  model  was  considered  adequate  to  facilitate 
comparison  of  autopilot  proposals  for  the  drone  boat. 

The  hull  shape  of  the  eight  metre  drone  boat  is  to  be  geometrically  similar  to  the 
six  metre  prototype.  Therefore  it  is  expected  that  the  manoeuvring  model 
developed  would  be  applicable  to  the  eight  metre  hull  after  dimensional  scaling  of 
the  coefficients  derived  for  the  prototype  boat. 


7.  Conclusion 

This  report  describes  a  simple  mathematical  model  of  the  prototype  RAN 
precursor  mine  sweeping  drone  boat.  The  model  was  based  upon  the  literature 
and  describes  a  boat  travelling  along  its  surge  axis  at  a  constant  speed  produced 
by  a  constant  thrust  from  a  single  outboard  motor.  The  model  predicts  the  yaw 
rate  of  the  boat  in  response  to  small  changes  of  the  angle  of  thrust. 


The  formulation  of  the  hydrodynamic  forces  acting  on  the  boat  is  based  on 
partial  derivatives  of  the  components  of  force  with  respect  to  the  dynamic 
variables.  The  manoeuvring  behaviour  of  the  vessel  is  described  by  these 
derivatives. 

A  sea  trial  was  held  at  Jervis  Bay  to  collect  manoeuvring  data  for  the  prototype 
drone  boat.  The  data  were  used  to  derive  coefficients  for  the  model.  These 
coefficients  are  considered  to  be  accurate  enough  to  provide  reliable  estimates  of 
drone  boat  performance  in  still  water  conditions. 
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10.  Nomenclature 


A  set  of  moving,  orthogonal  right  hand  Cartesian  coordinate 
axes  fixed  to  the  boat,  at  some  arbitrary  origin  O. 

Arbitrary  origin  fixed  to  the  boat. 

A  set  of  orthogonal  right  hand  Cartesian  coordinate  axes  fixed 
in  space  to  some  arbitrary  origin. 

Resultant  externally  applied  force  vector  acting  on  the  boat  at 
the  centre  of  gravity. 

Resultant  externally  applied  moment  vector  acting  on  the 
boat  about  the  origin  O. 

The  resultant  externally  applied  moment  vector  acting  on  the 
boat  about  the  centre  of  gravity. 

Linear  velocity  vector  of  the  boat  with  respect  to  the  centre  of 
gravity. 

Rotational  velocity  vector  of  the  boat  about  the  centre  of 
gravity. 

Position  vector  in  local  coordinates  of  the  boat  centre  of  gravity 
with  respect  to  the  local  origin  O. 

Rotational  velocity  vector  of  the  boat  about  the  local  origin  O. 

Linear  velocity  of  the  boat  at  the  local  origin  O. 

Moment  of  inertia  about  the  local  origin  O.  The  moment  of 
inertia  is  given  by  a  (3x3)  matrix. 

Boat's  moment  of  inertia  about  the  centre  of  gravity.  The 
moment  of  inertia  is  given  by  a  (3x3)  matrix. 

Time. 

Surge  Force  along  the  x-axis. 

Sway  Force  along  the  y-axis. 

Moment  about  the  z-axis. 

Steady  surge  velocity. 

Surge  velocity. 

Steady  surge  acceleration. 


Surge  acceleration. 


Steady  sway  velocity. 
Sway  velocity. 


V. 

v 

• 

Vo  Steady  sway  acceleration. 

• 

V  Sway  acceleration. 

Steady  yaw  rate, 
r  Yaw  rate. 


To 

T 

V 

m 

I„ 

x0 

8 

Yt 

Y„ 

Ye 

nt 

Y, 

N. 

▼ 

Y. 

▼ 

Nr 

Yr 

N. 


Steady  yaw  acceleration. 

Yaw  acceleration. 

The  instantaneous  boat  heading. 

Dry  mass  of  the  vehicle. 

Moment  of  inertia  about  the  z-axis. 

Position  of  the  boat's  centre  of  gravity  along  the  x-nxis. 
Outboard  motor  angle. 

Resolved  component  of  thrust  forces  in  the  sway  direction. 

Resolved  component  of  hydrodynamic  forces  in  the  sway 
direction. 

Other  external  forces  which  act  on  the  body  apart  from  thrust 
and  hydrodynamic  forces. 

Damping  moment  derivative  with  respect  to  sway. 

Damping  force  derivative  with  respect  to  sway. 

Added  moment  of  inertia  coefficient  with  respect  to  sway. 

Added  mass  coefficient  with  respect  to  sway. 

Damping  moment  derivative  with  respect  to  yaw. 

Damping  force  derivative  with  respect  to  yaw. 

Added  moment  of  inertia  coefficient  with  respect  to  yaw. 
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Added  mass  coefficient  with  respect  to  yaw. 

Exciting  force  coefficient  with  respect  to  the  outboard  angle. 
Exciting  moment  coefficient  with  respect  to  outboard  angle. 
Cross-track  error. 

Points  defining  the  position  of  the  baseline  track  in  (xe  ,ya) 
plane. 
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